Time for pulse traversal through slabs of dispersive and negative (e, /i) materials 
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The traversal times for an electromagnetic pulse traversing a slab of dispersive and dissipative 
material with negative dielectric permittivity (e) and magnetic permeability (fi) have been calculated 
by using the average flow of electromagnetic energy in the medium. The effects of bandwidth of 
the pulse and dissipation in the medium have been investigated. While both large bandwidth and 
large dissipation have similar effects in smoothening out the resonant features that appear due to 
Fabry-Perot resonances, large dissipation can result in very small or even negative traversal times 
near the resonant frequencies. We have also investigated the traversal times and Wigner delay times 
for obliquely incident pulses and evanescent pulses. The coupling to slab plasmon polariton modes 
in frequency ranges with negative e or /i is shown to result in large traversal times at the resonant 
conditions. We also find that the group velocity mainly contributes to the delay times for pulse 
propagating across a slab with n = — 1. We have checked that the traversal times are positive and 
subluminal for pulses with sufficiently large bandwidths. 

PACS numbers: 42.25.Bs, 42.65.-k, 78.20.Ci, 41.20.Jb 
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I. INTRODUCTION 

The time for traverse of light through a dispersive 
medium is interesting and important, both from a fun- 
damental viewpoint P, Q , and for technological applica- 
tions such as designing delay lines or systems for en- 
hanced non-linear applications. However, there are a 
variety of time scales, depending on the physical quan- 
tity being measured, that can be defined for this time 
of traverse [3, [1, d, Hj. A popular measure for the delay 
time of pulses is the Wigner delay time[f| t w = ^\ u =q, 
i.e, the frequency derivative of the phase of the output 
wave evaluated at the carrier frequency u>. The Wigner 
delay time, which is based on tracking a feature on the 
pulse moving at the group velocity (v g = can turn 
out to be superluminal or even negative, while at the 
same time well describing the motion of pulses with nar- 
row bandwidths over short distances [1, 0,111] ■ While there 
is no contradiction of Causality or of Special Relativity 
in these phenomena involving holomorphic pulses, the 
Wigner delay time becomes inaccurate for large pulse 
bandwidths or when there is a large deformation in the 
pulses. 

A measure that is based on the flow of electromagnetic 
energy for the time of traverse between two points n, and 
rf was proposed by Peatross et al[9(, which is given by, 



Eq. (2) can be rewritten exactly as [§] 



At = (t) r{ - (t) t 
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represents the arrival time of a pulse at a point r. This 
time scale is particularly suitable for pulses with large 
bandwidths as the relative contributions due to the prop- 
agation at the group velocity (group delay), and defor- 
mation of the pulse (reshaping delay) can be identified. 
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which is very useful for spectral calculations. This def- 
inition of the traversal time based on the motion of the 
centroid of the Poynting vector has been supported by 
experiments on dispersive mediafic|, and even in angu- 
larly dispersive systems [Tl[. It has also been shown that 
the arrival times for pulses measured through the rate 
of absorption in an ideal impedance matched detector 
is equivalent to the above arrival times [ijj]. We have 
also shown earlier that the definitions of the group delay 
times, and the reshaping delay times get interchanged for 
evanescent pulses [12]. 

Negative refractive index media (NRM) or left handed 
media (LHM) simultaneously have Re(e) < 0, and 
Re(/i) < at a given frequency [TI|, and have captured 
the imagination of the Physics community by their nu- 
merous counter-intuitive electromagnetic properties (See 
for a review of NRM). In isotropic NRM, the wave 
vectors k, and the Poynting vector S point in oppo- 
site directions. NRMs are also necessarily dispersive 
and dissipative in nature, and can additionally also sup- 
port surface (plasmons) states on their interfaces with 
positive media. Thus the study of the phase velocity 
{v p = x = ~)> the group velocity, and the energy flow in 
such media is interesting, and important. In fact in some 
metamaterials, one can have all possible combinations of 
(positive, and negative) phase, and group velocities 

In this paper, we study the times for pulse traversal 
through slabs of dissipative media with negative mate- 
rial parameters [e < 0, \x < 0) using the average energy 
flow in the media given by [9j. Pulse propagation in 
NRM has been principally studied with a focus on neg- 
ative refraction at interfaces [HI, [ItJ , and nonlinear 
effects [H, EH Hil- The Wigner delay time has been 
studied for pulses at normal incidence on slabs of NRM 
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in the limit of zero dissipation [2l|. We have earlier in- 
vestigated the traversal times based on the energy flow 
in infinitely extended NRM, and for normally incident 
pulses in semi- infinite NRM [l2j]. There the geometry 
was chosen such that the traversal times were affected by 
only the intrinsic dispersion of the medium parameters, 
and avoided the effects of other resonances such as slab 
resonances (by using infinite or semi-infinite media), and 
surface plasmon resonances (by normal incidence of the 
pulses). Here we will consider these effects of resonances 
on the traversal times for the transport of narrowband, 
and broadband pulses through slabs of causal NRM. We 
also show that finite levels of dissipation qualitatively 
change the nature of the traversal times near the reso- 
nances. Coupling to the slab plasmon polaritons will be 
shown to give rise to large delay times. We find that 
for pulses with large enough frequency bandwidth, the 
traversal times are positive, and subluminal. 
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FIG. 1: A dispersive slab having thickness Ar (Region 2) 
surrounded by vacuum (Regions 1, and 3). The black circle 
in Region 1 just outside the slab represents the source of ra- 
diation. The black screen in Region 3 just outside the slab 
boundary represents the detector. 



II. TRAVERSAL TIME ACROSS A DISPERSIVE 

SLAB 

We calculate the propagation times for electromagnetic 
pulses across a dispersive slab. The relative dielectric 
permittivity and the magnetic permeability of the slab 
are considered to be plasma like and Lorentz dispersive 
in nature, and given by, 



UJ 2 — U)\ + ioj^n 
UJ 2 — + iu^fn 



(4) 



(5) 



Here to Pj b,o=^fp,b,o with / p , /;,, and fo taken to be 12 
GHz, 6 GHz and 4GHz respectively. The slab behaves 
as a positive refractive index medium (PRM) or right 
handed medium (RHM) (e > 0, fi > 0) when u> > w p , a 
negative refractive index medium (NRM or LHM) (e < 0, 
fi < 0) within loq < uj < u>b, and as a barrier (e < 0, 
[i > 0) elsewhere. Here we note that the expressions 
of e and /i are similar to those given by [2l[ with the 
main difference arising due to the introduction of finite 
amounts of damping in the medium. 

For convenience, we take the source of radiation to be 
placed in vacuum just outside one boundary of the slab 
and the detector just outside the other one (Fig. 1). We 
have taken same medium (vacuum) on either sides of the 
slab (i.e. regions 1, and 3). So here E\ = /Xi = £3 = 
A*3 = 1- £2, and /j.2 have same forms respectively given 
by Eqs. (4), and (5). We take our initial pulse form to 
be Gaussian in time which is represented by, 

t 2 

E(r, = 0, t) = £E cxp[ -1 exp(-iwt). (6) 

Hence the Fourier transform of the pulse is given by, 
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where <D is the carrier frequency and r is the pulse du- 
ration. The magnetic field is simply obtained using the 
Maxwell's equations, 



H(r i ,w)=y— -= re 

2V2 



(8) 



where k z \ represents the wave vector in the first medium. 
We consider P-polarized light for normal incidence, ft 
should be noted that for normal incidence, the results 
are independent of the state of polarization of radia- 
tion whereas they are dependent in case of oblique in- 
cidence. Hence for the latter case, we will deal with 
both the S- and the P- polarizations separately. For the 
P-polarization, the magnetic field at the detector is re- 
lated to that at the source via the transmission coefficient 
across the slab. The final magnetic and the electric fields 
at the detector are given by 

Hfr^l^^-^re-^^TH, (9) 
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where k z ^ represents the wave vector in the third 
medium. 

Here T(w) represents the transmission coefficient 
across the slab which is given by, 



T(uj) = 4 



tt'e ik ^ Ar 

^ _ r /2 e 2ik z2 Ar ' 



(11) 



where Ar represents the slab thickness and k z i represents 
the wave vector inside the dispersive slab, t, t', and r 1 
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respectively represent the Fresnel coefficients of transmis- 
sion, and reflection by the slab interfaces and are given 
by, 



(RUM) 



(LHM) 



2^2 

£2 



k z2 

£2 



kz3 
S3 



kzl I kz2, ' 
£l £2 



S-2 



£3 



>^2 

■?2 



£3 



Here the unprimed, and primed coefficients stand respec- 
tively for the coefficients across the first and the sec- 
ond boundaries. For S-polarization, in the expressions 
of the Fresnel coefficients, the e's are simply replaced by 
fi's. Also the Fresnel coefficients relate the electric fields 
across the interface rather than the magnetic fields. Suf- 
fixes 1, 2, and 3 respectively represent the parameters 
at the source, slab, and the detector sides as described 
earlier. For convenience we later substitute equal ma- 
terial parameters on both (source and detector) sides of 
the slab. We calculate the delay times for different band- 
widths of the pulses. For this we use the same notation 
for the broad and the narrowband pulses as discussed in 
our earlier paper [l2j . The pulse has a broadband when 
wr=10 or less and it has a narrowband when ujt= 100 
or more. 

The Wigner delay time was calculated using Eq. (11), 
and is given by, 
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for P-polarization and 
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for S-polarization. Throughout our calculations (both for 
normal and oblique incidence), we have checked that the 
Wigner delay time yields the same result as the traversal 
time for narrowband pulses with the average energy flow 
method. 
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FIG. 2: (Color online) Scaled total delay time of a pulse with 
different bandwidths (wr) plotted as a function of frequency 
(f) across a dispersive slab with large thickness (Ar = 10cm). 
(a) Solid, and dashed lines show delay times respectively for 
a narrowband, and a broadband pulse across a nondissipative 
slab of RHM. (b) Same as (a), but across a slab of LHM. 
(c) Delay time for a narrowband pulse across a moderately 
dissipative (fy p — 0.01lo p , 7 m = O.OlWf,) slab (solid line), and 
a highly dissipative (7^ = 0.1w p , 7 m = O.luob) slab (dashed 
line) of RHM. (d) Same as (c), but across a slab of LHM. 
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A. Traversal times for normal incidence 

In this case, the parallel component of the wave- vector 
is zero and the pulse is normally incident on the slab. 
So there is no coupling with the slab plasmon polaritons. 
Since fc x = 0, the Maxwell's equations can be combined 
to give, 



(15) 



This is independent of whether the slab is of a RHM or 
a LHM. In Fig. 2, we plot the delay times scaled with 



FIG. 3: (Color online) (a) Scaled total delay time as a func- 
tion of the carrier frequency (f) for a narrowband pulse (solid 
line), and a broadband pulse (dashed line) across a highly 
dissipative (7 P = 0.1uj p , 7 m = O.luJb) slab of LHM, and large 
thickness (Ar = 10cm). (b) same as (a), but with a small 
slab thickness (Ar = 1cm). 



the free space propagation (— ) versus the carrier fre- 
quency for both broadband (Cut— 10) and narrowband 
(wt= 1000) pulses. For narrowband pulses, we refer to 
[21j | where the Wigner delay times were calculated for a 
nondissipative slab and it was shown that resonant fea- 
tures appear in the delay time behaviors due to presence 



of the poles of the transmission coefficient (Fabry-Perot 
resonances). We have taken (Ar = 10cm) as large thick- 
ness and (Ar = 1cm) as small thickness of the slab rel- 
ative to the wavelength (2.5cm) of the pulse at the elec- 
trical plasma frequency (f p ) . 

First we compute the results for the traversal times of 
broadband and narrowband pulses through a nondissipa- 
tive slab which is achieved by substituting 7 p = 7 m = 
in the expression of e and /i. Figs. 2(a) and 2(b) respec- 
tively show the traversal times for pulse propagation in- 
side slabs with positive and negative refractive indices. In 
both the figures, it can be observed that, the features due 
to the slab resonances get smoothened with an increase in 
the pulse bandwidth. So it is expected that for extremely 
broadband light, these features might completely disap- 
pear. Here we note that, the results for the narrowband 
pulses in Figs. 2(a), and 2(b), are exactly same as those 
for the Wigner delay times given in [2l[. Next we study 
the traversal times for narrowband pulses propagating 
through dissipative slabs of both RHM and LHM (Figs. 
2(c), and 2(d)). To include moderate levels of dissipation 
in the medium, we use 7 P = 0.01u p and 7 m = O.OloJb 
and for high levels of dissipation, we use 7 P = 0.1aj p and 
7 m = O.luib respectively in Eqs. (4) and (5). We see 
that when a small amount of dissipation is introduced 
in the medium, the time taken for transmission through 
the slab is less than that taken for the nondissipative 
case. With increased dissipation in the slab, one can also 
clearly observe that the slab resonant features just disap- 
pear. We have also studied the delay times for a highly 
dissipative slab of LHM for both narrow and broadband 
pulses (Figs. 3(a), and 3(b)). Fig. 3(a) shows the results 
for a slab with large thickness (10cm) and 3(b)shows the 
corresponding results for a slab with a small thickness 
(lcm). It can be clearly seen that the delay time is very 
small near the magnetic resonance frequency. Then it 
rapidly increases for large frequencies, and after passing 
through a peak, it gradually decreases. For narrowband 
pulses, the total delay time near the resonance frequency 
(u)q) even becomes largely negative. Even for broadband 
pulses, with small thickness of the slab, this negativity in 
the delay time appears near u>o although to smaller and 
smaller extent with increasing thickness of the slab. The 
anomalous dispersion of the refractive index of a medium 
with high amount of dissipation leads to small/negative 
delay times for broadband/narrowband pulses near the 
resonance frequency. It can be seen that, for narrow- 
band pulses, and large thickness of the slab, the delay 
time bends down near LOb where the value of \x becomes 
0. For the broadband pulse, the traversal time remains 
unaffected at u>b- 

Next we investigate the delay times within the fre- 
quency range (u>b < u> < uj p ) where the medium behaves 
as a plasma and as a consequence, most of the wave com- 
ponents are evanescent. We study the traversal times 
for both broadband and narrowband pulses. We find 
that for different slab thicknesses (Ar =5cm, Ar =10cm, 
Ar =20cm) involved, and regardless of the amount of dis- 
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FIG. 4: (Color online) (a) The scaled total delay time of 
a broadband pulse as a function of the carrier frequency 
(f) across a nondissipative slab which behaves as a plasma 
medium with different thicknesses shown by dashed line 
(Ar = 5cm), solid fine (Ar = 10cm), and dotted line 
(Ar = 20cm). (b) same as (a), but the pulse has a nar- 
row bandwidth, (c) same as (a), but the siab is highly dis- 
sipative (7 P = O.lojp, 7 m = O.lujb). (d) Fabry Perot reso- 
nance structures near resonance frequency in the delay time 
graph for a broadband pulse across a slab with large thickness 
(Ar=10cm). Dotted, solid, and the dashed lines indicated by 
numbers 1, 2, and 3 respectively show the features across 
slabs of nondissipative — 0, 7 m = 0), moderately dissi- 
pative (jp — O.OlLUp, 7 m = O.Olu;,), and highly dissipative 
(7p = O.lujp, 7 m = 0.1u>f,) media. 



sipation in the slab material, the traversal time remains 
the same for almost all frequencies except at the edges 
of the stop bands (Figs. 4(a), 4(b), and 4(c)). Such an 
independent behavior of the traversal time with distance 
as well as dissipation is the famous Hartman effect (22J. 
We also see Fabry-Perot resonance-like structures with 
large peaks appearing near the edges of the stop bands 
frequencies for broadband pulses with large thickness and 
nondissipative slab material (Fig. 4(d)). Such resonant 
features do not appear if the slab thickness is reduced, 
or interestingly, if the bandwidth of the pulse is small 
either. We can understand this surprising resonant effect 
for broadband pulses by noting that the plasma acts as a 
spectral filter. Consider a pulse with a large bandwidth 
and the carrier frequency within the stopband but near 
the lower edge. In transmission, most of the spectral 
components within stopband are almost completely at- 
tenuated while the spectral components in the passband 
at lower frequencies will get transmitted. If the band- 
width of these propagating spectral sub-components is 
small, then the behavior of the transmission as a func- 
tion of frequency will be modulated by Fabry-Perot reso- 
nances. Note that the impedance mismatch between the 
slab and vacuum will be much larger at u>b, thus empha- 
sizing the resonant effects. Such resonant effects disap- 
pear for small slab thickness because there the bulk of the 
spectral components within the stopband are not attenu- 
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FIG. 5: (Color online) Total delay time of an extremely narrowband pulse (lot — 5000) plotted versus scaled carrier frequency 
(f/f P ), and scaled parallel wave vector (q x ) for P-polarization across a moderately dissipative (~/ p = 0.01w p , j m = O.Olo;;,) slab 
with small thickness (Ar = -^-). (b) same as (a), but across a slab of highly dissipative (7^ = 0.1lu p , 7 m = OAlu^) material 
and large thickness (Ar = 2A P ). (c) Fabry Perot resonance structures in the delay time plotted versus scaled carrier frequency, 
and parallel wave vector across a moderately dissipative (-y p = 0.01lu p , 7™ = O.Olcjf,) slab with large thickness (Ar = 2A P ). (d) 
same as (c), but the slab is highly dissipative (7^ = 0.1u) p , 7™ = O.lwf,). 



ated enough for the Fabry-Perot resonances to dominate. 
Increasing dissipation also obviously destroys the contri- 
bution of these resonances to the phenomena. It should 
be noted that the transmitted pulses will be significantly 
stretched or deformed as their bandwidth will be sbstan- 
tially smaller. We have not studied here this aspect via 
the deviation (Ar) 2 of the delay time. 



B. Traversal times for oblique incidence 

Next, we study the traversal time of a pulse when it is 
obliquely incident on a slab, i.e., the parallel wave vector 
(k x ) is nonzero. For this case, the Maxwell equations 
give, 

k z = ~ k l = (16) 

where k x — q x ^. 

We consider the incident pulses with either S or P po- 
larizations. For convenience, we have scaled all the fre- 
quencies with respect to the plasma frequency. So here 
f p = 1, /o = 0.33, and /& = 0.5. The thicknesses of 
the slab used are Ar = -g- (small thickness), for which 
fcpAr = 0.4-7T, and Ar = 2X p (large thickness), for which 



fcpAr = 4ir. The thicknesses are relative to the wave- 
length (A p ) at the plasma frequency. This is important to 
note that, here for narrowband pulses, we use Cut — 5000 
(extremely narrowband case) and for broadband pulses, 
Cut — 1 (extremely broadband case). 



1. Propagating waves 

Here we discuss the pulses for which the wavevector is 
real, or in other words, the waves are propagating. First 
of all, we present the results for an extremely narrowband 
pulse (Cut — 5000) with moderate amounts of dissipation 
in the slab material. The delay time is plotted versus 
both the frequency, and the parallel wave vector (Figs. 
5(a), 5(b), 5(c), and 5(d)). From the graphs, it can be 
clearly observed that violent dispersion of the delay time 
occurs at the magnetic resonant frequency for both the 
polarizations. For P- polarization, a small dip in the 
delay time also occurs at the electric plasma frequency 
(f p ) for small thickness of the slab (Fig. 5(a)), and at 
both the fp, and fb for large thickness of the slab. For 
S- polarization, this additional small dispersion in the 
delay time curve apart from that at /o occurs only at the 
magnetic plasma frequency where as it becomes smooth 
at the electrical plasma frequency for all thicknesses of 




FIG. 6: (Color online) (a) Total delay time of an extremely broadband pulse (lot = 1) plotted versus scaled carrier frequency 
(f/f P ), and parallel wave vector (q x ) across a moderately dissipative (7^ = 0.01lo p , j m — 0.01a;;,) slab with small thickness 
(Ar — ^-), and P-polarization. (b) same as (a), but the slab is highly dissipative (7,, = 0.1lo p , 7 m = O.liot). (c) same as (b), 
with large slab thickness (Ar = 2A P ). (d) same as (a), but for S-polarization. 



the slab. This arises because the Fresnel coefficients for S- 
polarized waves do not depend on the dielectric constants 
(e) as strongly as the magnetic permeability (/1). When 
there is high amount of dissipation in the material of 
the slab, it is seen that, the scale of the delay time axis 
as well as the sharpness of the dispersion at the resonant 
frequency decrease simultaneously with the broadening of 
the peaks for both the polarizations. Here we should note 
that, the traversal times for narrowband pulses across 
slabs with large thickness and high amount of dissipation, 
give exactly same features with the same scales for both 
the S-, and P- polarizations. 

Then the delay time results are analysed for narrow- 
band pulses and large thickness of the slab in the fre- 
quency range which causes positive refractive index. Af- 
ter introducing moderate amount of dissipation in the 
medium, one can see Fabry-Perot resonance structures 
appearing beyond f p (Fig. 5(c)). Such resonant struc- 
tures (ripples) disappear with increase in dissipation in 
the medium (Fig. 5(d)). 

Next we study the traversal times for an extremely 
broaband pulse (lot — 1). It can be seen that the de- 
lay times are large at frequencies close to fo and large 
wave vectors (Figs. 6(a), 6(b), 6(c) for P-polarization), 
and (Fig. 6(d) for S-polarization). This happens for all 
thicknesses of the slab with any amount of dissipation in 
the material. In the delay time graph for S-polarization, 
a large peak can be seen at frequency close to fo and large 



wave vector. This occurs due to the large dispersion in fx 
at the magnetic resonance and the magnetic plasma fre- 
quencies. The delay time decreases with increase in dissi- 
pation in the medium. It can be observed that the graph 
becomes flat for small thickness of the slab, whereas for 
large slab thickness, the delay time increases rapidly at 
frequencies close to fo and large wave vectors. For large 
thickness of the slab with high amount of dissipation in 
the material, both the polarizations give similar features 
with equal scales (Fig. 6(c) for P-polarization). We 
should note that, for broadband pulses, the total time 
taken for traversal is positive and subluminal for all the 
parameters studied here. 



2. Evanescent waves 

Here we consider pulses for which most of the wave 
vectors are imaginary making the incident wave evanes- 
cent. This is achieved by making the second term in Eq. 
(16) under the square root large than the first term. We 
calculated both the Wigner delay time and the Energy 
delay time for such pulses most of whose components are 
evanescent. While the Wigner delay time can be easily 
calculated using the phase shifts, calculating the traver- 
sal time via the energy flow for evanescent pulses is a 
non-trivial problem. This is because the energy flow as- 
sociated with a single evanescent wave in vacuum is zero. 
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FIG. 7: A layered structure consisting of five slabs used to 
study the tunneling of pulses. Region 3 is the dispersive slab 
with thickness Ar with £3, and [13 respectively given by Eqs. 
(4), and (5). Regions 2, and 4 are air slabs (e.2 = A*2 = 1), and 
with thickness d. Regions 1, and 5 are semi-infinite slabs with 
£1 = £5 = 25, and y,\ = /X5 = I. The black circle in Region 
1 just outside the boundary of slab 2 represents the source, 
and the black screen in Region 5 outside slab 4 represents the 
detector. 



Thus the poynting vector for pulse with all spectral com- 
ponents having evanescent wave vectors is zero at the 
detection point in vacuum. Hence it is not possible to 
compute the energy traversal time for evanescent pulses 
if the slab is embedded in vacuum. One needs to couple 
the energy in these systems out to measure the pulses. 
This is related to the measurement of tunneling quantum 
particles whereby one needs to raise the particles with 
negative energy (evanescent waves) above the barrier be- 
fore detection. Similarly, we make an arrangement with 
layered slabs where the evanescent waves are outcoupled 
to propagating modes in high-index media. 

In this arrangement, we have taken two nondispersive 
slabs of different parameters kept symmetrically on ei- 
ther sides of the dispersive slab making a layered slab 
structure (Fig. 7). The first and fifth slabs have semi- 
infinite extent with large relative dielectric permittivity 
(e = 25) and relative magnetic permeability (/i = 1). The 
second and the 4th slabs consist of vacuum with e = 1, 
and \x = 1 and large slab thickness with corresponding 
kpAr equal to (47r). The 3rd or the middle dispersive 
slab has e, and \i respectively given by Eqs. (4), and (5) 
and small slab thickness with corresponding k p Ar equal 
to (0.47r). The source is present in the first medium just 
outside the boundary of the second slab and the detector 
is placed in the fifth medium just outside the boundary 
of the 4th slab. The value of q x in Eq. (16) is chosen in 
such a manner that the wavevector is real, making the 
pulse propagating in first and fifth slabs, and imaginary 
making the pulse evanescent in second and 4th slabs. 

First, we plot the Wigner delay time versus both the 
frequency and the wave vector in a moderately dissipative 
slab. In a rather uniform landscape of delay times, the 
resonant conditions for the slab surface plasmon polari- 
tons (SPPs) stand out in stark contrast where the mag- 
nitude of the delay times are comparatively very large. 
Thus the entire dispersion of the SPPs can be traced out 



(Fig. 8). There are two distinct plasmon modes cor- 
responding to the symmetric and antisymmetric modes 
whose frequencies tend to at large wave vectors. Simi- 
larly two modes also appear below the magnetic resonant 
frequency. For highly dissipative slabs also, such plas- 
mon modes are seen for evanescent waves, but with large 
broadening of the dispersion curves (Fig. 8(b)). The sur- 
face plasmon features are lost when the thickness of the 
slab is larger than X p . We have seen that the features 
of the plasmon modes for the S-polarization is different 
from the P-polarization as in that case, these modes of 
magnetic character occur between /o (magnetic resonant 
frequency), and fb (magnetic plasma frequency) when 
fi < 0, and superpose at 

Then we studied the traversal times for evanescent 
pulses having narrow bandwidth (lot — 5000) using the 
energy transport method with our new arrangement of 
the layered slab structure (Fig. 7). For this arrangement, 
we plotted the delay times for narrowband pulses versus 
both the frequency and the wavevector, and analysed the 
results for moderate and large dissipative slabs with both 
P-polarization (Figs. 8(c), and 8(d)) and S-polarization. 
We see that Figs. 8(c), and 8(d) look almost same as 
Figs. 8(a), and 8(b). Thus, the energy traversal times 
are also significantly affected at the Surface Plasmon Po- 
lariton frequencies. It is worth noting that the traversal 
times are large at the resonant conditions. 



3. Propagation through a slab with n = — 1 

Finally, we consider a slab having unit negative re- 
fractive index (n = —I) and surrounded by vacuum 
(n = +1). Negative refractive index of unit magnitude 
can be achieved at a single frequency for a nondissipa- 
tive slab. The properties of such a slab with (n = —I) 
are very interesting due to the possibility of designing 
a perfect lens [H, [23| By choosing f p = 1, / = 0.33, 

f b = 0.5, we get n = -1 at / = ^ {s=-l, /i=-I). 
With propagation inside the medium, the propagation 
distance increases by a factor of . Here due to perfect 
impedance matching, no multiple reflections take place. 
Using the expression for group velocity (v g — fe), the 
group delay along the direction of propagation is given 
by G d = -^s = -^(n + wl^). For the particular 

J u 7J„ cos r. cos 6 V- ouj ' 



frequency / 



the second term within the bracket in 



V q COS 9 

- Ll 

V2' 

the above expression gives a value equal to y-. We plot 
both the Wigner delay time, and the group delay time 
versus q x for propagating pulses for the particular fre- 
quency mode described above (Fig. 9). From the graph, 
it can be observed that the delay time gradually increases 
with q x until q x = 1 (where it becomes infinity). It can 
be seen that the graph feature of the group delay time 
is very similar to the Wigner delay time. Hence it is 
inferred that the group delay mainly contributes to the 
total delay occuring during the propagation of a pulse 
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FIG. 8: (Color online) (a) The Wigner delay times for pulse traversal with P-polarization as a function of frequency and parallel 
wave vector across a moderately dissipative (7^ = O.OltJp, 7m = 0.01a;i,) slab with small thickness (Ar = (b) Same as (a), 
but for a highly dissipative (y p — 0.1uj p , 7 m = O.lo;;,) slab. (c)The traversal times obtained by the energy transport method 
across the layered slab structure where the middle dispersive slab is moderately dissipative (7^ = O.OltJp, 7,„ = O.Qluib) in 
nature with small thickness (Ar = -§■"). (d) same as (c), but for a highly dissipative (7^ = O.lup, j m = O.lwj,) slab. The 
dispersion of the slab plasmon polariton modes of the slab stand out clearly and the resonant conditions for these modes are 
characterized by large Energy delay times. 



inside a slab with n = — 1. 



III. CONCLUSIONS 

In summary, we have demonstrated the manifold im- 
plications of dissipation on the traversal time of a pulse 
across a dispersive slab. Throughout our calculations, 
we have used average energy flow method [U to obtain 
the delay time of a pulse, and have checked numerically 
that the results obtained using the Wigner delay time 
method, and those using the average energy flow method 
are exactly same for very narrow bandwidth pulses. In 
our results we have shown that high amount of dissipa- 
tion in slab material, along with large pulse bandwidth 
smoothen out the resonant features. 

We have analysed the reason behind the occurrence of 
small/negative delay times near the magnetic resonant 
frequency which is a consequence of anomalous disper- 
sion of the refractive index of the slab medium. We have 
shown that large slab thickness, along with high material 
dissipation give rise to same features for both S and P 
polarizations in the case of oblique incidence. We have 
also shown that the group delay mainly contributes to 
the total delay across a slab with unit negative refraC- 



^S 15 



— Wigner delay time 
x Group delay time 




0.18 0.38 0.58 0.78 

% 



FIG. 9: (Color online)The Wigner delay time (solid line), and 
the Group delay time (cross symbol) plotted versus scaled 
parallel wave vector^) across a slab with small thickness 
(Ar = and unit negative refractive index (n = —1) sur- 
rounded by vacuum. 



tive index, and surrounded by vacuum. We have checked 
that the total time taken for a broadband pulse is always 
positive, and subluminal. 
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